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1. Introduction 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix ^ is a set of fundamental parameters in 
the Standard Model, which relates the mass eigenstates and the weak eigenstates in the electroweak 
theory. In the Wolfenstein parameterization it is given by 
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Because of unitarity, it contains only 4 independent parameters, (A,A,p,T]). 

To determine each CKM matrix element, one requires both theoretical and experimental in- 
puts. On the theoretical side, one needs to know relevant hadronic amplitudes, which often contain 
nonperturbative QCD effects. A major role of lattice QCD is to calculate such hadronic amplitudes 
nonperturbatively, from first principles. One can then extract the CKM matrix elements by com- 
bining lattice QCD results as the theoretical input with the experimental input such as branching 
fractions. 

To accurately determine each CKM matrix element from lattice QCD, one should use hadronic 
processes whose amplitude can be reliably calculated with existent technique. There are a set of 
such processes (amplitudes) — the so-called "gold-plated" processes (quantities) which contain 
at most one hadron in the initial and final states [Q]. These include the exclusive semileptonic 
decays and leptonic decays of B, D and K mesons, and neutral B — B and K — K mixings. For 
such processes, technique for lattice simulations are already well established, and thus reliable 
calculations are possible. 
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Figure 1 : Gold-plated processes for each CKM matrix element . The neutral K — K mixing (characterized 
by the CP-violating parameter Ek) is another gold-plated process, which gives a constraint on the phase of 
the CKM matrix (p , rj ). 

The gold-plated processes for each CKM matrix element are summarized in Fig. [T[ The mag- 
nitude of the CKM matrix element, e.g., \Vcd\ can be determined from either the semileptonic 
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decay D Tilv or leptonic decay D Iv, as explained below. \Vub\ and \Vtd\ can be respectively 
extracted from the semileptonic B — > ttZv decay and neutral B — B mixing, which give constraints 
on the phase of the CKM matrix (p, r]). The neutral ^ — ^ mixing gives another constraint. Taking 
these together, one can extract (p, T]) assuming that the Standard Model is correct. In this way one 
can, in principle, determine all 9 CKM matrix elements, i.e., all 4 Wolfenstein parameters using 
lattice QCD. 

The accuracy of the CKM matrix elements from lattice QCD is subject to several systematic 
uncertainties, however. The most serious one is the uncertainty from the "quenched" approxi- 
mation, in which effect of virtual quark loops (dynamical quarks) is neglected ("rif = 0"). This 
approximation has been adopted in the community for a long time, simply to reduce the computa- 
tional cost. Recent developments of computer resources and algorithms enable us to perform more 
realistic lattice calculations — "unquenched" simulations which include effect of light (up, down) 
dynamical quarks ("«/ = 2") or light and strange dynamical quarks ("«/ = 2+ 1"). For the current 
status of the unquenched simulations, see Ref. 

Most of gold-plated quantities listed above have been or are being calculated in unquenched 
lattice QCD. Given this situation, it would seem timely to present the CKM matrix elements from 
lattice QCD in a uniform fashion in one place. This paper gives a result for the whole CKM matrix - 
all 9 CKM matrix elements and all 4 Wolfenstein parameters - determined from lattice QCD using 
recent results for gold-plated quantities. For this purpose I exclude results from quenched QCD. 
However, I note that quenched calculations are still important and useful to study methodology 
and other systematic uncertainties. For recent reviews on the quenched calculations of gold-plated 
quantities, see Refs. ||, 0, |, |, |^]. 

Although there are several ways to discretize quarks on the lattice, up to now some of the 
gold-plated quantities have been calculated in unquenched QCD only using the staggered-type 
fermion. ^ This is because the staggered fermion is computationally much faster than other lattice 
fermion formalisms such as the Wilson-type fermion, domain wall fermion and overlap fermion. 
As a consequence, the results for the CKM matrix elements in this paper are often estimated from 
only one or two unquenched calculations. I hope that more unquenched results using other lattice 
fermion formalisms will appear in the near future, and leave future reviewers to make a more 
serious average of the CKM matrix. 

To present the result for the CKM matrix in a uniform fashion, I use lattice QCD results only 
as the theory input for nonperturbative QCD effects. It is of course desirable to include non-lattice 
theory inputs (such as ones for inclusive B decays) to improve the precision, but it is beyond the 
scope of this paper. For recent progress on non-lattice approaches for CKM phenomenology, see. 



e.g., Refs. |12,hJ, 14|. 



In the staggered fermion formalism, one needs to take the fourth-root of the fermion determinant to adjust the 
number of quark flavor, as a consequence of the fermion doubling. The validity of this procedure is not yet proven, so 
further study on this issue is necessary. For a review, see Ref. [O]. 
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The results for the magnitudes of the CKM matrix elements are 
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The results for the Wolfenstein parameters are 

A = 0.2249(17) , 
A = 0.77(7) , 
P = 0.16(7) , 
n = 0.37(4) , 



(1.2) 



(1.3) 
(1.4) 
(1.5) 
(1.6) 



where p = p(l - AV2) and fj =r]{\-X^/2). 

The rest of this paper is organized as follows. In Section ^, I review recent results for the 
semileptonic and letonic decays of B, D and K mesons, and determine the magnitude of the CKM 
matrix elements using them. 5 CKM matrix elements are directly determined using lattice results, 
whereas other elements (except for \Vtd\) are indirectly determined using CKM unitarity. In Sec- 
tion ^, I review recent results for the B° — B" and — mixing amplitudes, and extract the CKM 
phase (p, T]) by performing the unitarity triangle analysis. In Section ^, I give the conclusion. 

The first attempt of the full determination of the CKM matrix using lattice QCD was made in 



Ref. [15], where only results from semileptonic decays were used. The determination of the CKM 
matrix elements using lattice QCD has been discussed for a long time; see, e.g., Refs. [|, [l^ for 
recent reviews. A more detailed review on recent lattice calculations of K meson physics (such as 
the — mixing and K nlv decay) can be found in Ref. [16]. 



2. CKIM magnitude from lattice QCD 

I start this section with examples of the determination of the magnitude of the CKM matrix 
elements. The CKM magnitude is often determined from the semileptonic decays (such as D ^ nlv 
and K Tllv), and in some cases, leptonic decays (such as D ^ Zv and K Iv) also provide an 
independent determination with comparable precision. Below I take \Vcd \ as an example to explain 
how to extract the CKM magnitude. 

The branching fraction of the semileptonic decay D — > nlv is given by 
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(2.1) 



where q = pj) — pj^ is the momentum transfer, <7niax = {mo — nijc)^ and /+ is the form factor defined 
below. The relevant hadronic amplitude is conventionally parametrized as 
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Figure 2: D meson decays relevant to the CKM matrix element |Ve£;|. The left panel is the semileptonic 
decay D nlv, and the right is the leptonic decay D Iv. 



where is the vector current involving the heavy and light quarks. By combining the lattice 
calculation of f+{q^) with the experimental measurement of the branching fraction, one can extract 
\Vcd\- Similarly |Vc.v|, \'Vub\, and \'Vus\ can be respectively extracted from D Klv, B Dlv, 
B Klv and K nlv, as listed in Fig. |l]. 

The branching fraction of the leptonic decay D ^ Zv is given by 



Br(D^/v) = (known factor) x/o^ I 



(2.2) 



where is the D meson decay constant defined through 



(2.3) 



Here is the axial vector current involving the heavy and light quarks. Lattice results for /o can 
be used to determine \Vcd\ once the branching fraction is measured by experiment. The branching 
fraction of D ^ / v (Dv — > Zv) decay is being (will be) measured by CLEO-c [ 17, 18, 19|] and that of 
K ^ Iv has been measured by many groups pO|]. Thus they can be used to extract \Vcd\ (|Vci|) and 
|V,„ | respectively. On the other hand, B ^ Iv decay is still difficult to measure since the branching 
fraction is suppressed by a small factor |Vu^,p. 



2.1 \Vcd\ and \Vcs\ from D meson decays 

There are (at least) two reasons to study the D meson decays in lattice QCD. One is to extract 
the CKM matrix elements, as mentioned above. The other is that, taken the CKM matrix elements 
from elsewhere, one can test lattice QCD by comparing lattice results for the form factor y^^'^(^) 
and decay constant fo with corresponding experimental results. In particular, if lattice results 
agree with experiment for the D meson physics, one can have confidence in similar quantities for 
the B meson physics such as /f and fs, which are phenomenologically important and cannot be 
obtained by other means. 
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Figure 3: Light quark mass dependence of the D ~+ n form factor in Ref. pi]]. Symbols are data points, and 
lines are chiral fits using the staggered chiral perturbation theory (solid lines) and ones using a linear ansatz 
(dashed lines). The dashed vertical line indicates the physical ud quark mass. 



The Fermilab/MILC collaboration presented the first unquenched (nj = 2 + 1) calculation of 



form factors of D ^ nlv and D Klv decays [21, 22] using the MILC "coarse lattice" ensemble 



(a-i « 1.6 GeV) [El. We use an improved staggered fermion action [^J] for light {u,d,s) quarks 
and an improved Wilson fermion action [ |25| ] with the Fermilab interpretation [^] for the charm 
quark. To combine the staggered-type light quark (1 -component spinor) with the Wilson- type 
heavy quark (4-component spinor), we convert the staggered quark propagator into the naive quark 
(4-component spinor) propagator according to 

^{x)\x{x)x{y)My) = {¥{x)¥iy)), (2-4) 

where D.{x) = i^'fyi^i^ ^ proposed in Refs. [|^, ^]. Although the naive quark propaga- 
tor describes 16 (=2^) equivalent fermions, known as the "doubling problem", only the phys- 
ical mode contributes to the low energy physics and the remaining 15 (= 16 — 1) modes de- 
couple when it is combined with the Wilson-type (or any doubler-free) heavy quark. This can 
be understood by noting that the heavy-light meson mass for each mode is roughly given by 
~ {m„ + m^, niu + {pic + 2r/a), lUu + {mc + 4r/a), • • •} with r being the Wilson parameter 
for the heavy quark action. By taking the asymptotic limit in the time direction {t —>■ oo) for heavy- 
light 2-point and 3-point functions, the state with the lowest mass (physical mode) can be isolated. 
This method has been successfully applied to the heavy-light meson physics by the HPQCD col- 
laboration and Fermilab/MILC collaboration.^ 



Note that this method may be applied to the K meson physics as well, with the staggered u,d quarks and the 
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Figure 4: Form factors of D ^ nlv and D—^Klv decays innf = 2+1 lattice QCD by the Fermilab/MILC 
collaboration [|T]]. The errors are statistical only. Recent experimental results by the BES collaboration 
and by the FOCUS collaboration S] are also shown. 



Since the staggered fermion is fast and free from the exceptional configurations, one can per- 
form simulations with the light quark mass mj as low as m jS, in contrast to previous calculations 
using other fermions. This situation is shown in Fig. |[ Physical results at m/ = m^d = {mu + nid) /2 
are obtained from chiral fits using the staggered chiral perturbation theory (S;^PT) [ p9| , ^ ], but a 
simple linear fit gives consistent results. This suggests that the results at the physical light quark 
mass are insensitive to the fit ansatz and the error from the chiral extrapolation (m/ mud) is under 
control. The availability of data at small light quark masses is crucial. On the other hand, since our 
calculation is done at a single lattice spacing, the error from the lattice discretization effects is still 
large, giving about 10% total systematic uncertainties for the form factors. 

Our final results are shown in Fig. ^ together with recent experimental results [31 , 32, 3^, p4| ]. 



Wilson-type s quark. 
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The results are obtained using the parameterization of Becirevic and Kaidalov (BK) [35], 



/+(^')=(r-2)(T3^' /o(^') = (2.5) 
where ^ = /m^t. We obtain [ pT[ 22] 

fO-^ = 0.64(3) (6), a^^'^ = 0.44(4) (7), jS^^'^ = 1.41(6)(13), (2.6) 
F^-^ = 0.73(3)(7), a^^^ = 0.50(4) (7), j8^^^ = 1.31(7)(13), (2.7) 

where the first errors are statistical and the second systematic. The results agree with experiment 



for both the normalization at 17^ = [Bjl, B3] and the 17"^ -dependence [B3l B4p. This may indicate 



2 



reliability of lattice results for the similar quantities for B physics, the B — > nlv form factor. By 
integrating out f+{q^) in terms of and using the experimental measured branching fraction [^], 
we obtain 

|V,rf|semi-iep = 0.239(10) (24) (20) , (2.8) 
iK.lsemi-iep = 0.969(39) (94) (24) , (2.9) 

where the first errors are statistical, the second systematic, and the third are the experimental errors 
from the branching fractions. 

The Fermilab/MILC collaboration has also finalized the nf = 2 + 1 calculation of the D and 



Ds meson decay constants [36, 37, 38]. The employed lattice actions are the same as the ones for 
the semileptonic calculations. We performed partially quenched simulations, where the valence 
light quark mass can be different from the dynamical light quark mass, at three values of the 
lattice spacing. The chiral extrapolation is done using the S;^PT formula [|3^, and the final results 
are obtained by taking 6a^ after the chiral extrapolation (upper figure of Fig. §), where da^ 
denotes constants in the S^PT which parametrizes lattice discretization effects from the staggered 
fermion. After 5a^ — > 0, the lattice spacing dependence is small (lower figure of Fig. Our final 
results are fo = 201 (03) (17) MeV and/o^ = 249(03) (16) MeV, where the first errors are statistical 
and the second systematic. The two largest sources of the systematic error are the discretization 
effects and the chiral extrapolation (for fo). 

The CP-PACS collaboration reported a new unquenched (nf = 2) calculation of the Dj^.) decay 
constant [ pO| , pTj ]. For the light quarks, they used an 0{a) improved Wilson action. For the charm 
quark, they used a relativistic on-shell improved action |^, ^] which is similar to one in Ref. [26] 



but derived from a different point of view in Ref. [|42|]. Their simulated light quark mass ranges 
> ms/2, and a linear chiral extrapolation was made. They studied the lattice spacing dependence 
using three lattice spacings as in the Fermilab/MILC calculation, and performed a continuum ex- 
trapolation combining two data sets (one is from the temporal axial vector current and the other is 
from the spatial current) as shown in Fig. |6[ Their preliminary results are fo = 202(12)(l25) MeV 
and fo^ = 238(1 l)(l27) MeV, where the systematic errors (second parentheses) are dominated by 
uncertainties from the continuum extrapolation. 

The chirally improved actions were applied to the charm quark physics in two recent quenched 
(uf = 0) studies. Ref. [ p4[ ] used the overlap fermion action, whereas Ref. [45] used a variant of 



the domain-wall fermion action. The use of the chirally improved actions has an advantage that 
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Figure 5: D and Ds meson decay constants in «/ = 2 + 1 lattice QCD by the Fermilab/MILC col- 
laboration [|3^. The upper figure shows the light quark mass dependence and chiral extrapolation for 
/Dj Jmo, I {fD^/mo). The lower figure shows the lattice spacing dependence of 0^ = fosJniD,. 



they are free from 0{a,mQa) discretization errors without any tuning of parameters, and so may be 
useful for the simulations involving the charm quark where niQa < 1 . 

Turning to the experimental result for the leptonic decay, the CLEO-c collaboration updated 
their measurement of D ^ /xv branching fraction |[r|, |l^]. Assuming that \Vcd\ = | | = 0.224(3), 
they obtain the experimental result f^^^^^^ = 223(17) (03) MeV, which is in agreement with the 
recent lattice results. Their precision, 0(10%), is similar to that for the lattice results. The agree- 
ment may indicate reliability of lattice calculations of the heavy-light decay constants. They will 
further improve the precision of //> and report the result for fo^ in the future. 

One may conversely use the CLEO-c result for the D ixv branching fraction to determine 
\Vcd I . If it is combined with the nf = 2 + 1 lattice result for fjj by the Fermilab/MILC collaboration, 
one obtains 

|Krf|iept = 0.250(22)(21), (2.10) 

where the first error is one from the lattice calculation and the second is from the experimental un- 
certainty. \Vcd\ from the leptonic decay is consistent with the value obtained from the semileptonic 
decay Eq. (^^, and the size of uncertainties is similar to each other. 

It is interesting to consider the ratio of the leptonic and semileptonic decay branching fractions 
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Figure 6: Lattice spacing dependence and continuum extrapolation (a — > 0) for fo, in nf — 2 lattice QCD 
by the CP-PACS collaboration Filled circles are results from the temporal axial vector current and 

open circles are from the spatial one. 



because the CKM matrix element \Vcd\ cancels in the ratio. Writing 



R 



I Br(D Iv) 
Br(D nlv) ' 



(2.11) 



the R is proportional to fo/U dq^iff^'^iq^))^]^^^. Since one can directly compare lattice results 
with experiment for this quantity, it provides a good test of lattice QCD. The experimental result 
(dominated by the CLEO-c measurements) is [|T^ 



/?exp = 0.25(2), 



and the unquenched lattice result using the Fermilab/MILC calculation of and f'_^_ 

/?lat = 0.21(3). 
They are in reasonable agreement with each other. 



IS 



(2.12) 



(2.13) 



As for theLattice'05 value of the CKM matrix elements, I take a weighted average of Eqs. (2.8) 



and ( |2.10| ) for \Vcd\, and simply quote Eq. ( p!9| ) for |Vci|; 



l^rfkatOS 
|Vci|Lat05 



0.245(22) , 
0.97(10) , 



(2.14) 
(2.15) 



where the errors are the combined uncertainties from theory (from lattice QCD) and experiment. 
These are consistent with the values quoted in Particle Data Group (PDG) [20]; |Vc£/|pdg = 0.224(12) 
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and |Vci|pDG = 0.996(13). Note that the PDG values above are obtained from neither the semilep- 



tonic nor the leptonic decays, and thus Eqs. ( [2. 14| ) and ( [2.15| ) provide an independent determination 
of \Vcd\ and |V„|. 



2.2 B meson decays 
2.2.1 B nlv decay 

The semileptonic decay B —>■ nlv can be used to determine \Vuh\, which is one of the most 
important CKM parameters to constrain the unitarity triangle. Since the pion momentum available 
in lattice calculations is limited up to around 1 GeV, only the higher ^ -region iq^ > 16GeV2) can 
be simulated for the B nlv decay. This is in contrast to the D — > nlv decay, for which we can 
cover all -region {0 < < g^ax)- Since the experimental accuracy of the branching fraction is 
better for the lower 17^ -region than the higher -region, it is desirable to extend the lattice result to 
the lower -region. Below I first summarize recent lattice results for the B —>■ nlv form factors, 
and extract |V„fo| using the results for the higher -region. I then discuss recent attempts to extend 
the results to the lower -region. 

There are two nf = 2+1 unquenched calculations of the B —>■ nlv form factors using the MILC 
coarse lattice ensembles with improved staggered light quarks. One is obtained with the NRQCD 
heavy quark 



16p by the HPQCD collaboration [47, 48], and the other is with the Fermilab heavy 
quark by the Fermilab/MILC collaboration [22, i9]. Their results are shown in Fig. ^. The total 
systematic uncertainties are 11% for the form factors in both calculations. The largest quoted error 
comes from the perturbative matching between the lattice and continuum theory for the HPQCD 
result, and from the discretization effect for the Fermilab/MILC result. Two unquenched results 
agree with each other within errors. These are also consistent with previous quenched calcula- 



tions [|50|, |51|, Qj, |53|, q4[] ; at present it is difficult to estimate the effect of the quenching quantita- 
tively because the size of other systematic uncertainties, 0(10%), is the same as the expected size 
of the quenching error. 



By integrating f+{q ) over 16 GeV < q < q^^ and using an average [55] of the partial 



branching fractions Br(16 GeV^ < r < measured by CLEO [56], Belle [57, 58, 59] and 



BABAR [6C, 51,^] collaborations, the CKM matrix element is obtained as 



IK^Impqcd 

iKfolpNAL/MILC 



4.47 (22)(49)(30)xl0"^ 
3.78 (30)(42)(25)xl0"^ 



(2.16) 
(2.17) 



where the first errors are statistical, the second systematic, and the third are the experimental errors. 
As for the Lattice 2005 value, I take a simple average of the two preliminary results, obtaining 



|V„i|Lat05 = 4.13(62)xl0 



-3 



(2.18) 



nlv form factors. Both HPQCD and Fermi- 



The total uncertainty for \Vuh \ is 15%. 

Let us discuss the ^^-dependence of the B 
lab/MILC collaborations use the BK parameterization Eq. (2.5) to interpolate and extrapolate the 
results in q^. To estimate the uncertainty from the BK parameterization, the Fermilab/MILC col- 
laboration also made a fit using a polynomial ansatz in q^. The difference between the two methods 
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Figure 7: Form factors of the B Ttlv decay inn f = 2 + 1 lattice QCD by the HPQCD collaboration @ Q 
and by the Fermilab/MILC collaboration [p2[ B9p. Lines are ones from fits with the BK parameterization. 



for \Vuh\'=< 



2^2^-l/2 



is 4% with q\ 



16 GeV^. This difference is included in the 



systematic error in Eq. ( [2.17 ). With ^^j^^ = 0, however, the difference amounts to be 11%, which 
is a significant effect. This is because a long extrapolation is required from lattice data points 
(16 GeV^ < < ^max) ^ = 0. For a more precise determination of |Ki|> it is necessary to 
reduce the uncertainty for the lower -region. 

One solution is to combine the lattice results for the higher -region with non-lattice results 
for the lower -region. Ref. | |63| ] combined the recent unquenched lattice results with the QCD 
dispersion relation and |y„fo|/+(0) from an analysis of the B^nn decay based on the Soft Collinear 
Effective Theory (SCET). The QCD dispersion relation is not a model but an analyticity bound. 
Ref. [54] made a similar study using /+(0) from the light-cone sum rule instead of one from 
the SCET. Ref. [65] combined quenched lattice results with the QCD dispersion relation and the 
experimental measured -dependence. In each case, the uncertainty for \Vub\ can be reduced by 
Ri5% with the additional information on the form factors. 

Another solution is a direct lattice simulation at lower q^ using the moving NRQCD (mN- 
RQCD) [B6l |63, BSl]. The mNRQCD is a generalized version of non-relativistic QCD in the B 
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Figure S: B ^ nlv form factors using the unquenched lattice results, QCD dispersion relation and 
|^Mfo|/+(0) from an analysis based on the soft colUnear effective theory 



meson moving frame (u = Pb/Mb / 0). The action is given by 

i^'mNRQCD = ( iD, + /(v • D) + - ^^-^ + • • • V (2.19) 

where = 7(1, v) and 7^^ = \/l — v^. Setting v = gives the usual NRQCD action. The mN- 
RQCD allows the B nlv calculation at lower with smaller p;t; for v ^ 0.75, = can be 
achieved with Pn = I GeV []69[], where the size of lattice discretization effects from non-zero pion 
momentum is modest. 

Previous lattice simulations with the mNRQCD for a large v suffered from large statistical 



errors []66|], but Ref. [ |70| ] showed that it can be reduced by using a special smearing function so that 
the simulation with v 0.7-0.8 is feasible. Ref. |^] made a quenched study to test the mNRQCD. 
They calculated the decay constant of a heavy-heavy meson /hh at non-zero momentum p using 
both the NRQCD (v = 0) action and mNRQCD (v / 0) action. As shown in Fig. |, /hh with the 
NRQCD depends on p, indicating that the discretization effect from non-zero pion momentum is 
not under control. On the other hand, the result with the mNRQCD is constant in p as it should 
be, and the statistical errors are reasonably small. The result with the mNRQCD looks encourag- 
ing, and so it may be worth studying the B —>■ nlv from factors at lower q^ using the mNRQCD 
formalism. 



2.2.2 B Dlv and B D*lv decays 

The form factors of B ^ Dlv and B ^ D*lv decays can be calculated more accurately than 
those of heavy-to-light decays, due to the approximate symmetry between the initial and final states. 
The branching fraction of the B ^ D/v decay is given by 

Br{B Dlv) = \Vcb\^ f dw F{wf x (known factor), (2.20) 
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Figure 9: Momentum-dependence of the ratio of heavy-heavy meson decay constants ///// (p) / ///// (0) from 
NRQCD (crosses) and mNRQCD (squares) 



where w = vb-vd with vg = PB/niB and vd = Po/niD- The form factor of the B ^ Dlv decay at zero 



recoil limit, Ffi_>z)(w = 1), can be precisely determined by considering the double ratio [73, 74]; 

CWoB(j)C^i/oZ)(f) {D\Vo\B){B\Vo\D) 



CWoD(j)Cfi^oB(f) {D\Vo\D){B\Vo\B) 



(2.21) 



where €^^"^[1] and {D\Vo\B) are respectively the B ^ D three-point function and amplitude. 

The Fermilab/MILC collaboration has a preliminary ?iy = 2 + 1 unquenched result for the 



B — > Dlv form factor using the MILC configurations [22]. The light quark mass dependence of the 
form factor is mild and a linear chiral extrapolation was made, as shown in Fig. [T^. The unquenched 
result is F^^^(l) = 1.074(18)(16), which is consistent with the quenched result [73]. Combined 



with an average of experimental results for |Vfi|F(l) [pSQ, we obtain 

\Vct\L.m = 3.91(09)(34)x 10-2, (2.22) 

where the error from the lattice calculation (first error) is much smaller than the experimental one 
(second). 

For a more precise determination of |Vcfo|, the B D*lv decay should be used because the 
experimental uncertainty is smaller. As for the lattice calculation, the chiral extrapolation is crucial 
because a singularity should appear for the form factors at mps = mo* — mj) which is close to the 
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Figure 10: Light quark mass dependence of the unquenched B Dlv form factor at zero recoil [J22[]. The 
quenched result []73[] is also shown around at m/ = 



physical pion mass point. Ref. [75] calculated the B — > D*Zv form factors in the staggered chiral 
perturbation theory (S;^PT), and found that the singularity (seen in the continuum ;tPT) disappears 
due to lattice discretization effects. The unquenched calculation of the S — > Z)*/v form factors is 
underway [[75|]. 



2.3 K meson decays 

The semileptonic K nlv decay is traditionally used to determine \Vus\- The form factor of 
K nlv at <7^ = can be precisely calculated using the double ratio method as in the B Dlv 
case. The first precise lattice calculation of the K nlv form factor has been done in the quenched 



approximation using improved Wilson quarks |76]. The quenched result is 



/f-'^(0) = 0.960(5)(7) inf = 0), 



(2.23) 



which agrees with an earlier estimate using a quark model [|77[]. The uncertainty is smaller than 
1% due to the approximate symmetry between initial and final states, as for the B Dlv form 
factor. Since then, three unquenched calculations have been started. One is the = 2 calculation 



using improved Wilson quarks by the JLQCD collaboration [|78|], obtaining the preliminary value 
/+ "^(0) = 0.952(6). Another «/ = 2 calculation is underway using domain-wall quarks by the 
RBC collaboration [fTOj]; their preliminary result is /f^'^(O) = 0.955(12). A preliminary = 
2 + 1 calculation on the MILC configurations is done using a combination of improved Wilson 
and staggered quarks by the Fermilab collaboration [15], getting f^^^{0) = 0.962(6)(9). These 



calculations rely on chiral perturbation theory to guide the chiral extrapolation (m/ ntud) using 
data at m/ > ms/2. It may be interesting to perform a direct simulation at smaller m/ using staggered 
fermions or chirally-improved fermions. For technical details of the calculations, see Ref. [p^]. 
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Taking a simple average of three preliminary unquenched results, I obtain 

/f-'^(0) = 0.956(12) {nf>2). 
Combining this with an experimental result for |Vi,.s|/^^^(0) [ pO| ] gives 

|V«.|semi-lep = 0.2264 (28)(12), 
where the first error is from lattice calculations and the second from experiment. 



(2.24) 



(2.25) 



\Vus\ may be determined from the leptonic decay K Iv []81[]. The lep tonic decay constant /k 
has been precisely calculated in = 2 + 1 lattice QCD using improved staggered quarks by the 



MILC collaboration [p3[]. Their updated result is p2| ] 

fK/f. = i.m{3)tll). 

Using this and the experimental result for Bi{K Iv) /Br(7r Iv) 

|Ks|iep = 0.2242(+i}), 



(2.26) 



one obtains 



(2.27) 



which is consistent with that from the semileptonic decay Eq. (2.25). 

For the Lattice 2005 value, I take a weighted average of Eqs. (2.25) and ( 2.27] ), obtaining 



|V,v|Lat05= 0.2249(17). 



(2.28) 



There is a 2a disagreement between Eq. ( [2.28| ) and the PDG average \Vi,s\ = 0.2200(26) pO|]. The 
reason is as follows. The PDG used the K ^ nlv decay for iKsl; the form factor is taken from 



Ref. [77 1 which is consistent with Eq. (2.24), but the K — > kIv branching fraction is from earlier 



experimental measurements which disagree by « 2a with the recent one used for Eq. (2.25) 



2.4 Other CKM magnitudes from CKM unitarity 

Having determined the 5 CKM matrix elements, one can check a unitarity condition of the 



CKM matrix using results from lattice QCD alone. From Eqs. (2.14), (2.15) and (2.22), one gets 

{\Vcd\^ + \Vcs\^ + \Vcb\^y^^ = 1.00(10), (2.29) 

which is consistent with CKM unitarity. 

Conversely one may use CKM unitarity to determine other CKM matrix elements as follows; 



|K^|Lat05 = {l-\Vus\^-\Vub\^y^^ 

= 0.9744(4), 

|Vf,,|Lat05 = (1-|K/,|'-|K*|^)^/^ 

= 0.9992(1), 

|V,,|Lat05 = \V:,Vub + V*,Vct\ I \M ^ \y*cych\ I \M 

= 3.79(53)xl0"2_ 



(2.30) 
(2.31) 
(2.32) 
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One can also determine some of the Wolfenstein parameters from Eqs. (2.28), (2.22) and (2.18). 
One gets 

ALatos = IK.I = 0.2249(17), (2.33) 
ALatos = iK-fol/A^ = 0.77(7), (2.34) 
(p' + '7')[{?05 = \Vuh\/{AX') = 0.47(8). (2.35) 

The determination of remaining CKM parameters, (p,?]) and \Vtd\, will be discussed in the next 
section. 

3. CKM phase from lattice QCD 

In this section, I extract CKM phase parameters (p , T] ) and the magnitude \Vtd \ using recent lat- 



tice results. One constraint on (p, t]) has already been obtained in the previous section, Eq. ( |2.35| ). 
Two more constraints may be obtained from the mixing of neutral B and K mesons. 

The neutral B mixing is characterized by the mass difference of 5" and B^ mesons, which is 
given by 

AMb,, = (known factor) x /IBb^ \V,lVt,\^ {q = d,s) (3.1) 
The nonperturbative QCD effects are contained in /J Bb^, where /b^^ is the B^ meson decay constant 



defined in an analogous way to Eq. (2.3) and Bb^ is the Bg meson bag parameter defined through 



{B'>\{bq)v-A{bq)v-A\B'') = -mB,,BBjl {q = d,s), (3.2) 

where {bq)v-A is the V — A current involving b and q quarks. By combining the lattice calcu- 
lation of IIjBbj with the precisely measured value of AMbj, one can extract ^ \Vtd\ = 
AAV(l-P)^ + '7^. which gives a constraint on (p , Tj ) . Since some uncertainties for lattice results 

cancel in the ratio J , it is desirable to consider the ratio of and B^ meson mass differences. 



Up to now, however, only a lower limit is known for AMb^. The CDF and DO experiments are 
expected to measure AMb/, once it is available, it will provide a better constraint on (p, T]). 
The neutral K mixing is characterized by the CP-violating parameter Ek, given by 

|%| = Bjf T][(l-p)ci + C2] (3.4) 

where ci and C2 are numerical constants, and Bk is defined through 

{K%d)v-Aisd)v-A\K^) = l^niKBKfl (3.5) 

with fK being the K meson decay constant. Combining the lattice result for Bk and the experimental 
result for |£/f| gives another constraint on {p,T\). 
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As for other constrains on (p,T]), recent experiments by B factories enable the measurement 
of all 3 angles of the CKM unitarity triangle (a,jS,7) = (arg (^^) ,arg ,arg 
In particular, the accuracy of j8 is impressive. See, for example, Refs. [83, 84, |l^] for recent 
experimental measurements of the CKM angles. 



The constraints, Eqs. ( 2.35 ), (3.1 ) and (3^), as well as other constraints, over-determine (p, Tj). 
If good precision can be achieved for each sector, one can test the Standard Model by seeing 
whether or not there is inconsistency between them. Below I review recent lattice results for /b, 
Bb and Bk, and then extract (p, T]) using them. 

3.1 B meson mixing (/bBb) 

3.2 /b 

Let us recall that /b is similar to fo and lattice results for the latter agree with the experimental 
result by CLEO-c within f«10% uncertainties, as seen in previous section. This may indicate 
reliability of lattice results for /b . 

This year the HPQCD collaboration finalized their rif = 2 + I calculation of /b using the 
improved staggered quarks and the NRQCD heavy quark on the MILC configurations ||8f 



88]. They performed simulations at the light quark masses in the range m,v/8 < mi < nis/l. The 
chiral extrapolations are made using various fit forms including ones with the staggered ;tPT, the 
continuum ;^PT and a simple linear ansatz. They reported that the various fits agree with each 
other within 3% after the extrapolation, suggesting that the results at the physical light quark mass 
(m/ = niud) are insensitive to details of the chiral fit. 



The Hf = l + \ result by the HPQCD collaboration is shown in Fig. |11| together with arif = 2 
result by the JLQCD collaboration [|9|]. Although the JLQCD result is consistent with the HPQCD 
result at m/ ~ rris/l, a linear chiral extrapolation using only JLQCD data (with m/ > ms/2) clearly 
deviates from the HPQCD result at m/ = m,,^/. I believe that this is evidence that data at m/ < ms/2 is 



required for high-precision lattice calculations. The result by the HPQCD collaboration is [85, 88] 



fB, = 216(09) (19) (07) MeV, (3.6) 
/b, = 260(07) (26) (09) MeV, (3.7) 

where the first errors are statistical ones and uncertainties from the chiral extrapolation, the second 
are ones from the 1-loop perturbative matching for the current renormalization, and the third are 
other uncertainties. The 2-loop or nonperturbative matchings will be required for an accuracy better 
than 5%. The uncertainties from the matching (and some others) cancel in the ratio of B^ and Bg 
decay constants, giving 

/b7/b. = 1-20(3)(1). (3.8) 

This is a 3% determination, thanks to the small error from the chiral extrapolation. Once Amg^ is 
measured, this will significantly reduce the uncertainty for iVfil/lV/^/l. 

Comparing the HPQCD results with the latest averages of unquenched results («y > 2) in 
Ref. [pH], one sees a reasonable agreement for the individual /^^ and fs^, and good agreement for 



the ratio /bJ fsj, as shown in Fig. 12. The good agreement for the ratio is probably due to the 
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Figure 11: The light quark mass dependence of unquenched fss^/f^J ifBdy/f^Bd) by the HPQCD collabo- 
ration (nj = 2 + I, black symbols) [ |8^ and by the JLQCD collaboration (nj — 2, red symbols) [89|. Black 
(red) lines are the chiral extrapolations using the HPQCD (JLQCD) data. The dashed vertical line indicates 
the physical light quark mass (m/ = m^d). 



cancellation of some systematic uncertainties, as mentioned above. I also note that the averaged 



value in Ref. | |90| ] is estimated by including the effect of the chiral logarithm from ;^PT to guide 
the chiral behavior for m/ < nis/l. 

3.2.1 Bb 

Turning to Bb, there is no new or updated result in unquenched QCD this year.^ The nf = 2 
calculation by the JLQCD collaboration using an improved Wilson light quark and the NRQCD 



heavy quark [89] is still only the result from unquenched QCD. Their result is 



= 0.836(27)(+^^), 
BbJBb, = 1.017(16)(+^^). 

Combining this and the HPQCD result for gives 

= 244(26) MeV 



(3.9) 
(3.10) 



(3.11) 



■'while this paper was being completed, a new unquenched {nj = 2) calculation using domain- wall light quarks and 
an improved static heavy quark has been reported [^]. The results for BbJBb^ and fs, / fsj \/ Bg^ I^B^ are consistent 
with but larger than Eqs (3.1G) and (3.12). 
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Figure 12: Comparison of the B meson decay constants in nj = 2+1 QCD by the HPQCD collabora- 
tion [pSp and the unquenched world average in Ref. 



1.210(+^^: 



(3.12) 



where Bb is the renormalization group invariant bag parameter. These values should be compared 
with the previous average, e.g., Ref. [@ quoted fB^x/BBj = 214(38) and /bJ/bj^/BbJBbj = 



1.23(6). Equation ( 3.11 ) together with the experimental value of Ambj leads to 

|V,</|Lat05 =7.40(79)xlO-^ 



(3.13) 



which is consistent but smaller than the PDG value [^], |V,rf|pDG = 8.3(1.6) x 10 ^. Equa- 



tion (3.12) and the forthcoming measurement of Amb^ will give a 3^% determination of |Vf.s|/|V,rf|. 

The bag parameters are also studied in quenched («/ = 0) QCD using the overlap light quark 
action [^2|], and in perturbation theory using the twisted mass light quark action [ l93| , |94| ]. Using 
these actions has an advantage that operator mixings do not occur (or can be removed), which may 
lead a more precise calculation of the bag parameters in the future. For previous quenched results, 
see, e.g., Ref. [0]. 

3.3 K meson mixing (Bk) 



Three new studies of Bk inrif = 2 + \ unquenched QCD are reported this year [95, 96, p7|]. 



In particular, the HPQCD collaboration presented a preliminary value using improved staggered 
quarks on the MILC configurations [^], 



Bf^(2GeV) =0.630(18)(130)(34), 



(3.14) 



where the first error is statistical, the second is from the 1-loop matching which is again the largest 
error, and the third is other uncertainties. Ref. [^] also used an improved staggered fermion on the 
same configurations. On the other hand, Ref. [ p7| ] used the domain wall fermion for both valence 
and dynamical quarks. At present the latter two groups do not quote the physical value of Bk. 
A preliminary study in nf = 2 unquenched QCD using the unimproved Wilson fermion is also 



reported in Ref. |98]. For details of recent Bk calculations, see Ref. [16]. 
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Figure 13: Comparison of recent unquenched results for Bk IP5L y% IOC] and an average of quenched 



results |90|. The nf ~2 result by the UKQCD collaboration is obtained with an improved Wilson fermion 
action and one by the RBC collaboration is obtained with a domain wall fermion action [100|. The 
error for the result by the RBC collaboration is statistical only. 



A comparison of recent unquenched results [95, 99, IOC] and an average of quenched re 



suits [|90p is shown in Fig. 13. The unquenched results are consistent with the quenched average, 
but uncertainties are much larger. Given this, I do not take the average for the unquenched Bk, and 



simply use Eq. (3.14) as a representative of unquenched results in the following unitarity triangle 
analysis. 

3.4 Unitarity triangle analysis 

Let us now analyze the unitarity triangle using recent lattice results to extract (p,f]). As 



mentioned before, I use unquenched lattice results only as the theory input. Here I adopt Eq. ( |2.35D 
for (p2 + 772)1/2 from the B nlv form factor, Eq. for fB,\l^,, and Eq. for Bk- 

As the experimental inputs, I use AM^^ and |£^|. The B —>■ nlv branching fraction BT{q^ > 
16 GeV^) is also used to obtain Eq. ( 2.35 ). Using the lower limit of AMbJ AMb^ together with 
Eq. ( 3.12 ) for /bJ/bj^/BbJBb^ does not affect the result of this analysis for the reason given 
below. I also use the experimental result for sin(2j3) (=0.726(37) [^5|]) from B {cc)K^*^ decays 



to see its impact on the (p , T] ) determination. 

The unitarity triangle using lattice results together with the sin(2j8) constraint is shown in 



Fig. [TJ. The shaded regions indicate 1 a error bands. A difference between this and previous ones 
(e.g., one in Ref. [pO|]) is that the position of the AMb^ bound has moved to the right with a smaller 
uncertainty. Consequently, the bound from the lower limit of AMbJAMbj does not change the 
result for (p,T7). 

From the overlapped region in (p , f] ) plane, I obtain 



0.23(12), 
0.41(07) 



(3.15) 
(3.16) 
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Figure 14: Unitarity triangle using recent unquenched lattice results for constraints from |Vj,t| (red), AMb 
(green) and Ek (blue). The sin(2j3) constraint (yellow) is also shown. The shaded regions indicate la error 
bands. 



without the sin(2j3) constraint. Including the sin(2j3) constraint significantly improves the preci- 
sion, giving 

PLatos = 0.16(7), (3.17) 
r]Lat05 = 0.37(4), (3.18) 

which are consistent with the PDG values pO| ] p = 0.20(9) and fj = 0.33(5), and the accuracy is 
comparable, as shown in Fig. |l^. 

4. Conclusion 

This paper presents a full determination of the CKM matrix using recent lattice results for 
gold-plated quantities. To extract the CKM matrix elements in a uniform fashion, results from 
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Figure 15: Comparison of (p,J7) using recent unquenched lattice results (without and with the sin(2j3) 
constraint) and that quoted in Ref. 



unquenched lattice QCD are exclusively used as the theory input for nonperturbative QCD effects. 
The results for the CKM matrix elements are Eqs. ^J^, ( p^ , ( ^ ), ( g^ , ( |230| ), 

(p3Tl), (E32|) and (Ol). The results for the Wolfenstein parameters are Eqs. (|233|), (p3|), (Ol) 



and ( |3.18| ). These are summarized in one place, Eqs. (|1.2|)-( |1.6| ). 

At present, many unquenched results are obtained with improved staggered fermion actions. 
On the other hand, fewer unquenched results are obtained with other fermion formalisms (such as 
the Wilson-type fermion, domain wall fermion and overlap fermion), especially for the rif = 2 + 1 
case. Consequently, the results for the CKM matrix elements presented here are often estimated 
from only one or two unquenched calculations. I expect that more unquenched results using other 
lattice fermions will appear in the near future, leaving future reviewers to make a more serious 
average of the CKM matrix. 

The unquenched results for D physics (such as the D —>■ nlv form factor /f and the D 
meson decay constant fo) are in agreement with recent experimental results. This may increase 
confidence in lattice results for similar quantities for B physics (such as the B — > nlv form factor 
and /b). 

The typical accuracy of most of gold-plated quantities is 0(10%), being dominated by un- 
certainties from the lattice discretization effects, perturbative matching, and the chiral extrapola- 
tion. The last uncertainty especially applies to lattice fermions other than the staggered fermion. 
To achieve an accuracy better than 5%, simulations at smaller lattice spacings and smaller quark 
masses and higher-order matchings will be required. 



013/23 



Full determination of the CKM matrix using recent results from lattice QCD 



Masataka Okamoto 




Figure 16: Expected unitarity triangle with 5% accuracies for the B ^ nlv form factor and Bk and a 3% 
accuracy for fBjfBj \/BbJBbj- The AMb, measurement is also assumed here. 



A better accuracy (3% or less) is obtained for the B — > Dlv and K Ttlv form factors due to 
the approximate symmetry between initial and final states, and for the decay constant ratio fsj fs^ 
due to the cancellation of systematic errors. The latter will lead to a more precise constraint on 
(p, T]) once AMfi^ is measured by experiment. 

Assuming the AMg^ measurement and 5% (or better) accuracies for lattice results, the unitarity 
triangle will be something like Fig. 16. I hope that this will be realistic in next 5 years so that we 
can more precisely test the standard model using lattice QCD and be ready for new physics. 
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